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Bi-modular adaptive CDMA receiver 

The present invention concerns in general terms a reception method and device 
operating in code distribution multiple access (CDMA) mode. More precisely, when 
the propagation of the signal between the transmitter of a user and an array of 
reception antennae takes place on multiple paths, the present invention makes it 
possible to estimate the signal transmitted by the transmitter from the antenna signals. 

In a fixed or mobile telecommunication system, the propagation of the signal 
from a transmitter to a receiver can take place by multiple paths, each path being 
characterised by a delay, a direction of arrival, an attenuation and a phase rotation. If 
the telecommunication system enables a plurality of transmitters to establish a link 
with the same receiver, as is the case in mobile cellular telecommunication systems, 
the signals from the different mobile terminals intended for the base station are 
transmitted by different channels. Each transmission channel comprises a plurality of 
paths as far as the base station, each path being characterised by a delay, a direction of 
arrival, an attenuation and a phase rotation. 

Certain systems, such as the ones used in a direct sequence code distribution 
multiple access (DS-CDMA) technique must deal with high interference. The 
different users transmit at the same time in the same frequency band and the signals 
transmitted necessarily interfere with each other. The separation of the different users 
is effected by multiplying each symbol of the user by a spreading sequence peculiar to 
it, also for this reason referred to as the signature of the user, the different signatures 
ideally being chosen so as to be orthogonal. The spreading sequence frequency (chip 
rate) being greater than the frequency of the symbols, the signal transmitted by each 
user is distributed (or spread) in the frequency space. On reception, the signal of a 
user is separated by virtue of a filtering adapted to the corresponding signature. This 
filtering is also referred to as "despreading". The ratio between the band occupied by 
the spread signal and the band occupied by the information signal is referred to as the 
spread factor. 

The signatures employed must have good correlation properties, namely a very- 
pronounced auto-correlation peak and low inter-correlation values. 



The first of these two characteristics allows synchronisation of the sequence 
received. It is very useful when the transmission channel of a user includes several 
propagation paths. This is because each path can then be isolated by virtue of a 
filtering adapted to the signature and to the delay of the path. Advantage can be taken 
of the diversity of propagation within the channel in order to increase the signal to 
noise ratio on reception. To do this, a bank of adapted filters is used, separating the 
different paths, and the outputs thereof are combined. The most widespread 
combination is MRC {Maximum Ratio Combining), which consists of multiplying the 
signal output from each adapted filter by the conjugate of the complex multiplicative 
coefficient introduced by the chaimel on the path concerned. The resulting filtering 
operation is a filtering adapted to the equivalent filter of the channel. Because of its 
structure, the receiver thus formed is referred to as a rake receiver. Naturally, perfect 
separation of the paths takes place only if the auto-correlation is a Dirac. Ln practice, 
however, the separation is not complete and leaves multipath interference, which is 
also referred to as self noise. Fig. 1 depicts schematically a DS-CDMA system with 
K users. The data of a user k are spread in fi-equency by the corresponding signature 
in the module 100k before being transmitted over a chaimel 11 Ok having P paths. On 
reception, for a given user k, the signals being propagated according to the different 
paths i=l..P of the channel are separated by adapted filters 120k,i...l20k,P (only the 
battery of filters of the user k has been depicted) before being weighted by a set of 
complex coefficients ci./. The signals thus weighted are added (140k) and the 
resulting sum at the output of the rake receiver is subsequently detected in order to 
provide an estimation of the data of the user k. In the case of a downlink (links from a 
base station to a mobile terminal) the channels 1 to .ST are identical whilst they differ 
in the uplink (Hnks from mobile terminals to the base station). The first case can, 
from this point of view, be considered to be a particular case of the second. 

The second characteristic set out above guarantees a low level of interference 
between two distinct users. Nevertheless, there also, in practice, the inter-correlation 
between two signatures is rarely zero. This is notably the case in a so-called dazzle 
situation {Near-far effect) where a high-power signal received from a user interferes 
with the reception of a low-power signal coming from another user. Moreover, when 
the number of users is high, close to the spread factor, the sum of the interference of 
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the different users, low if taken in isolation, can have very disturbing effects for 
detection. 

More elaborate techniques of eliminating interference are then often necessary 
to mitigate the orthogonality defects of the sequences. To combat the interference 
5 between the channels of the different users (so-called multi-user interference), several 
methods have been proposed. A review of these will be found in the article by Simon 
Moshavi entitled "Multi-user detection for DS-CDMA communications" which 
appeared in IEEE Communications Magazine, Oct. 1996, pages 124-136. Amongst 
the existing multi-user techniques, the techniques of subtractive elimination 

10 (Subtractive Interference Cancellation) have good performance for reasonable 
complexity in use. These techniques apply well to the simple situation where the 
transmission channel of a user includes a single path. In this case, the filter modelling 
the channel can be limited to the multiplication by a complex coefficient. When the 
channels are multipath, the situation is on the other hand much more complicated 

15 since it is necessary to eliminate both the multipath interference and the multi-user 
interference. An iterative detector with subtractive elimination of the multi-user 
interference when there are multiple paths was proposed in an article by M.C. Reed et 
al. entitled "Iterative Multiuser detection using antenna arrays and FEC on multipath 
channels" published in IEEE Journal on Selected Areas in Communications, Vol. 17, 

20 N° 12, Dec. 1999, pages 2082-2089. Each iteration of the detector comprises a 
filtering adapted to each propagation path of the transmission channel of the user, a 
channel formation and a combination of the rake type. These three steps can be 
considered to define the treatment of a rake receiver modified in that it includes a 
formation of channels according to the directions of arrival of the different paths. 

25 Such a reception device is illustrated in Fig. 2. The device comprises an array of 

antennae (200i) ...(2000- Each antenna (200^) is connected to a battery of adapted 
filters (210j,k), i=l..P, each filter (210i) being adapted to a path i of the user k. The 
output signals relating to the same path i=l..P are directed to channel fomiers 
(220i,k,..,220p,k). Each channel former (220i,k) receives from a direcfion of arrival 

30 estimator (230i,k) the estimated value 6 of the angle of arrival 6^,^ of the path / of the 
transmission channel k, which enables it to point its reception beam in the 
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corresponding direction of arrival. In the case of an array of uniform linear antennae, 
the channel formation (220i,k) consists of effecting the calculation: 

yi,k=^xeAk.exp{-J^e,i,k) where ^ej,k=27td/A.cos0i.k.(£-l) (1) 

and where xi,i,k and yi.k axe respectively the inputs and output of the channel 
5 former (220;,^). 

0 i_k is the angle of arrival defined with respect to the axis of the linear antenna 
array, d is the distance between two antennae in the array and X the length of the 
wave transmitted. 

The outputs of the P channel formers are multiplied at (240i,k) by complex 
10 coefficients C/,jt before being added by an adder (250) in order to supply an output of 
the MRC (Maximum Ratio Combining) type. The complex coefficients c,- * are the 
conjugates of estimations of the complex multiplicative coefficients of the different 
paths i of the user K that is to say ci,k=ou,kJsyip{-jvi,k) where oajc and are 
respectively the estimated coefficient of attenuation and the estimated phase rotation 
15 of the path / of the user k. The result of the combination is therefore: 

Zk =^yuk. di,k.exp(-j Vi,k) (2) 

This reception device makes it possible to effect both a channel formation in the 
directions of the different paths and a combination of the MRC type as in a rake filter. 

The method proposed presupposes however that the coefficients of attenuation, 
20 the phase rotations and the directions of arrival of all the paths of all the users are 
determined. 



The aim of the invention is to propose a rake receiver including a channel 
formation but not necessitating the explicit calculation of the coefficients of 
attenuation, phase rotations and directions of arrival of the different paths. 
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This aim is achieved by the object of the invention, defined by a method of 
receiving a signal transmitted by a transmitter and arriving at an array of antennae 
after being propagated along a plurality of paths, comprising a filtering step 
decomposing each antenna signal into separate signals issuing from the different 
5 paths, a channel formation step forming the signals of paths from the said separated 
signals by means of a first set of complex coefficients, a combination step Hnearly 
combining the said signals of paths by means of a second set of complex coefficients 
in order to supply a combined signal, the method forming a plurality of first error 
signals between a reference value of the signal transmitted and the said signals of 
10 paths, and a second error signal is formed between the said reference value and the 
said combined signal, the first and second sets of complex coefficients being adapted 
so as to respectively minimise the root mean squares of the first signals and of the 
second error signal. 

According to a first embodiment, the signal transmitted comprising modulated 
15 data in the form of symbols, the reference value is a pilot symbol. 

According to a second embodiment, the signal transmitted comprising 
modulated data in the form of symbols, the method also includes a step of estimating 
the symbol transmitted from the combined signal, the said reference value then being 
an estimated value of the transmitted symbol. 

20 According to a third embodiment, the signal transmitted comprising modulated 

data in the form of symbols, the method also includes a step of estimating the symbol 
transmitted from the combined signal, a step comprising a demodulation of the 
estimated symbol into estimated data, a deinterleaving and a charmel decoding of the 
said data, a step comprising a channel recoding of the decoded data followed by an 

25 interleaving and a modulation of the said data in order to supply a re-estimated value 
of the transmitted symbol, the reference value then being chosen so as to be equal to 
the said re-estimated value. 

Advantageously, the first and second sets of complex coefficients are adapted 
according to an algorithm of the LMS or RLS type. 
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The object of the invention is also defined by a device for receiving a signal 
transmitted by a transmitter, the device comprising an array of antennae and 
comprising means adapted to implementing the steps of the method disclosed above. 

The characteristics of the invention mentioned above, as well as others, will 
5 emerge more clearly from a reading of the description given in relation to the 
accompanying figures, amongst which: 

Fig. 1 depicts schematically a DS-CDMA telecommunication system; 

Fig. 2 depicts schematically a rake receiver including a channel formation; 

Fig. 3 depicts schematically a receiver according to one embodiment of the 
10 invention. 

The basic idea of the invention is to dispense with an estimation of the 
directions of arrival, coefficients of attenuation and phase rotations of the different 
paths by initialising the complex coefficients of the channel formation and of the 
combination of paths at predetermined values and applying an algorithm for adapting 
15 these coefficients from the measurement of a signal of error between a reference 
signal and a signal observed after the recombination step. 

Let the situation be assumed where a user k transmits symbols Sk(t) over a 
transmission channel k consisting of a plurality of propagation paths i=l ..P. Each path 
i of the channel k arrives at an array of antennae at an angle of incidence 6i,k(t) . The 
20 array of antennae is preferably a linear array consisting of L antennae regularly spaced 
apart by a distance d, but other types of array, notably circular arrays, can also be 
envisaged. 

Let xtxk be the signal received by the antenna i , £=(l,. . .,L) corresponding to 
the path of index i of the user of index k\ 



25 



xtu(t)=[aiMt).&^V{J^'-'')\ ■ [e.-)i.v(j27vdlXcose,,k{t).{£-\ ))] .sk{t)+nuMt) 0) 
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The first term between square brackets represents the attenuation and phase 
rotation undergone by the signal on the path / of the channel k, the second term 
between square brackets represents the phase shift due to the difference in operation 
between the antenna i and the reference antenna 1 and nt,i,k(t) represents the white 
5 additive Gaussian noise at the antenna i , including the thermal noise, the multipath 
interference and the multi-user interference. The signal transmitted sk(t) then 
undergoes a double phase rotation resulting from the channel effect {vi,k(t)) on the one 
hand and the differences in operation between the antennae in the array on the other 
hand ( <pe.i,k(t)=27rd/lcosei,k(t).{i-\ ) ) . 

10 Fig. 3 illustrates a rake receiver including a channel formation with adaptation 

of the complex coefficients, as proposed by the invention. 

The receiver depicted is dedicated to the reception of the symbols transmitted 
over a channel k. The system has L antennae represented schematically by the block 
(300k). Each of the Z antenna signals is then correlated by a battery of P adapted 

15 filters where P is the number of propagation paths of the channel, each filter being 
adapted to a path i. The set of L.P adapted filters is represented diagrammatically by 
the block (310k), the L output signals, jcf,,-.* , i=\..L, of the filter adapted to the same 
path i being directed towards a channel former associated with this path. The P 
channel formers are depicted in the form of a single block (320k). The P output 

20 signals, yuk , of these channel formers are then multiplied by weighting coefficients 
and added in the block (340k) in order to supply a combination signal zk . The 
combination signal is submitted to a decision unit (360k) supplying, in soft or hard 
form, the symbol sk transmitted by the user k. Soft form means here the result of a 
weighting of symbols by their respective a posteriori probabiUties. The estimated 

25 symbol sk is then demodulated in order to supply one or more estimated data dk . 
These data can then be subjected to an operation of deinterleaving and channel 
decoding (370k) with soft or hard values (provided, naturally, that the data had been 
interleaved and coded by a channel coding corresponding to the level of the user). 
The channel former unit (320k) performs the following operations: 
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(4) 



where be,i,k are complex weighting coefficients. Operation (4) can then be 
written in an equivalent manner in the following vector form: 



yi,k(t)=hi,k(t) ^ ■%,k{t)=(h,i,k(t)bi.i.k(t} • -bLuit)). 



5 where be.i.k(t)=ae,i,k.ex.p(-j'¥(-'''d with i = (X....L) 

Secondly, the different paths of the signal of the user of index k are combined 
by the module (340k). The combination signal z* is written: 

zk(t)=Ck(t) ^ .yk(t)={a,k(t)cu(t} ■ <:p.k(0} 

where ci,k(t)=pi,k.cx.pi-C'rk) with / = 0 P) 

10 The invention proposes to adapt the coefficients bi.i,k so as to minimise, in root 

mean square value, the first error signals output fi-om the channel formation module 
(320) and to adapt the coefficients a,k so as to minimise, in root mean square value, a 
second error signal output firom the path combination module (340). 

The first error signals are expressed by: ^i,k(t)=q'i,k(t)-yi,k(t) (6) 

15 The second error signal is expressed by: ^k(t)=q"k(t)-Zk(t) (7) 

where q'i,k{t) and q"k{t) are reference values output from (320) and (340) 
respectively. The reference values q\kit) and q"k{t) are chosen so as to be equal: 

q\,k{t)=q"k(f)=qk{t) where qk(t) is a common reference value which does not depend on 
the path. 



X2.i,k(t} 



(4') 



(5) 
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To adapt the coefficients bi,i,k and Ci,k, it is possible to use, as the common 
reference value, as indicated symbolically by the switch (390k), either directly pilot 
symbols pk(t) or estimated symbols Skit) output from the decision module (360^), or 
the symbols s'k{t) obtained after channel decoding and receding. In the latter case, 
5 the estimated s3niibols Sk{f) are first of all demodulated and the data dk obtained are 
deinterleaved, decoded by channel decoding (370ic) and then the data thus obtained are 
receded by channel coding, reinterleaved and remodulated (380k) in order to supply 
new estimated symbols, s\{{) , having a better reliability than the initially estimated 
symbols Sk{t). Advantageously, the decisions on the symbols taken at (360k) and 

10 those on the data taken at (370k) have soft values, the symbols s'kif) then being 
obtained by the weighting of symbols in the constellation by their respective a 
posteriori probabihties. It should be noted that receding is not necessary if the 
encoded data aie recovered directly at the output from the decoding. For example, if 
an algorithm of the SOVA (Soft Output Viterbi Algorithm) type is used for decoding, 

1 5 the data in coded form can be recovered directly along the maximum likelihood path. 

The adaptations of the coefficients be,i,k and a,k can be effected with the pilot 
symbols only or be pursued with the estimated symbols, before or after channel 
decoding, or be effected solely by virtue of the latter in the case of complete absence 
of pilot symbols. 

20 

The coefficients hxk and a* are initialised at the start of the process of 
adaptation to the following values be,i.k{0')=S(£-£o) and a.k(0)=c for each path i where 
S is the Dirac symbol, £o is an antenna number and c a complex coefficient (for 
example equal to 1). The starting configuration is thus an omni-directional 
25 ' configuration (without channel formation) and multipath. 

Alternatively, the coefficients fo.a can be initialised by means of an estimator 
external to the device according to the invention, supplying estimated values ^a(O) of 
the angles of arrival and estimated values Vi,i(0) of the phase rotations of the different 
paths at time 0. More precisely, the initial values b(,i,k{0) are chosen such that: 
30 £)f,a(0)=exp(-y(i/a(0)+2;rc?/;i.cos4,A(0).(^-l))). The coefficients ca can also be 
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initialised by means of an estimator external to the device according to the invention, 
supplying estimated values ai,k(0) of the attenuations undergone by the signal along 
paths of the transmission channel at time 0. Then a,k(0)=di,k(0) is chosen. 

Whatever the case, during the adaptation, the coefficients are modified so as to 
5 minimise the mean values and . Because of the feedback of the error 

signal ^i,k , the values ifre,i,k converge progressively towards the values Vi,k+<pe,i,k and 
consequently compensate for the phase differences due to the charmel and to the 
difference in operation between anteimae. As a result the chaimel formers are no 
longer content, as in the configuration illustrated in Fig. 2, to form beams in the 
10 directions of the different paths but also, at least partially, compensate for the phase 
rotations caused by the different paths. The module (340k) then effects a combination 
on the signals of paths pre-compensated in phase issuing from the channel formers 
and corrects the phase difference remainders. 



1 5 square error. Many algorithms have been proposed for this purpose in the literature 
and a description of the main ones can be found in the jovunal article by L.C. Godara 
entitled "Application of antenna arrays to mobile communications. Part II: Beam- 
forming and direction-of-arrival considerations" published in Proceedings of the 
IEEE, Vol. 85, N° 8, August 1997, pages 1 195-1245. 

20 For example, it is possible to use an algorithm based on the stochastic gradient 



The coefficient adaptation algorithm is based on a minimisation of the root mean 



method. 



More precisely, if 




the coefficients are adapted as follows: 
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b/.K?+l)-ba(0-y"'.*-VA*(0 (8) and Ck{t+l)=Ckit)-^\,VJ'k(i) (9) 
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where and /^"^ are positive real coefficients defining the coefficient 
adaptation pitches. 

If the so-called LMS (Least Mean Square) algorithm is used in its complex form, 
the gradients VA* and VJ"a are estimated at time t by -£'a(0-x:,(0 and -£"k(t).yi{i) 
5 respectively, whence, taking account of (8) and (9): 

ba(/+l)=b/,KO+/^'u.£'a(0.x:,(0 (8') 

a(t+l)=Ck(t)+jLJ\.^'ic(t).rXt) (9') 

where the sjrmbol * designates the conjugation operator. 

To ensure convergence of the adaptation, the coefficients and fi"'k are 
10 chosen so as to be less than l/rr(Rf.t) and l/Tr(Rk) respectively, where Ri,k and Rk 
are the correlation matrices of the respective vectors Xf,i(0 and yk(t) and Tr(M) is the 
trace of the matrix M. 

Advantageously, coefficients ju'ik and //"jt can be chosen with low values after 
the algorithm has begun to converge, in order to increase the degree of precision of the 
15 convergence. Equally, it is possible to choose coefficients ju\k and //"^t with larger 
values after the coefficient has commenced converging, in order to increase the ability 
to continue the algorithm. It is also possible to use a standardised LMS algorithm, in 
which the adaptation pitch ) is divided by the norm of the input vector 

(x;,a(0 and yk(t) respectively) in equations (8') and (9'). 

20 As an alternative to the LMS algorithm, the adaptation can be effected by means 

of the RLS (Recursive Least Square) algorithm as described, for example, in the 
article cited above or by means of an algorithm derived from the latter, preferably the 
Fast-RLS algorithm as described in the article by D. Falconer et al. entitled 
"Application of Fast Kalman estimation of adaptive equalization" published in IEEE 

25 Transactions on Communications, Vol. COM-5, Oct. 1978 , pages 1439-1446. 
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The RLS algorithm is based on an inversion by recurrence of the correlation 
matrix ( Ri^k, Rk)- The complex coefficients are then adapted in the following 
manner: 

Ck{t+l)=ck(t)+€i,k(t).ki(t) (10) where kkit)=yk(t) ^.Pk(t) (11) 

5 where kt is an updating vector referred to as the Kalman gain and Fk(t)=R-^^(t) 

is the inverse of the correlation matrix calculated by the recurrence formulae: 

PA(0=PA(/-l)-P*(?-l).y,(Okjt(0 (12) 



with P*(0)=lp where Ip is the unit matrix 
and the Kalman gain ki(/) is calculated by: 



n(t) 

k,(0= ^ (13) 



where rk(t) is the linear form defined by n(t)=Pk{t-l).y^(t) (14) 

The adaptation of ba(0 is effected of course by equations similar to (10), (11), 
(12), (13), (14). 

15 Advantageously, a forgetting factor will be included in equations (12) and (14) 

so that only the contributions of the most recent input vectors will be taken into 
account in the estimation of the correlation matrix. 

According to one embodiment of the invention, the adaptation of hi,k(t) will be 
effected by means of the standardised LMS algorithm and the adaptation of Ck(t) will 
20 be effected by means of the Fast-RLS algorithm. 
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Although the invention has been disclosed in the context of a DS-CDMA 
system, it can apply to other access modes provided that the channels of the different 
users and the different paths taken for a given channel can be separated on reception. 
Where the signals received are in broad band, it is important to provide an alignment 
5 timewise of the signals which underwent different propagation times. If the 
separation of the different users and paths is effected by means of batteries of filters at 
the antennae, the outputs of the filters relating to the same path must be delayed so as 
to be aligned timewise. Likewise, prior to the recombination, an alignment in time of 
the different paths must be provided for. The alignment in time may be effected by 
1 0 FIFO buffers with mobile output pointers adjusted finely by interpolation. 

In addition, although the receiver of Fig. 3 has been depicted in the form of 
modules, it is however clear to a person skilled in the art that the operations 
represented by the different modules in Fig. 3 can be performed by a single processor 
progranomed for this purpose or by a plurality of dedicated processors. 

15 



